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This  report  details  an  investigation  into  the  development  of  in  situ  methane 
as  a  tracer  for  oxygen  uptake  in  plunge  pools  at  low-head  hydraulic  structures. 
The  research  was  conducted  by  the  St  Anthony  Falls  Hydraulic  Laboratory 
(SAFHL),  Department  of  Civil  and  Mineral  Engineering,  University  of  Minne¬ 
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Hydraulics  Laboratory  (HL),  WES.  Mr.  Steven  C.  Wilhelms,  RWQB,  was  the 
WES  point  of  contact.  The  work  was  conducted  under  the  general  direction  of 
Messrs.  Frank  A.  Herrmann,  Jr.,  Director,  HL,  and  Glenn  A.  Pickering,  Chief, 
Hydraulic  Structures  Division,  HL.  Dr.  John  W.  Keeley  was  Director,  EL. 
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availability  of  the  gas  chromatograph  and  to  Mr.  John  R.  Thene  for  advice  on 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  meters 

feet 

0.3048 

meters 

gallons  (U.S.  liquid) 

3.785412 

liters 

inches 

2.54 

centimeters 

1  Introduction 


Hydraulic  structures  have  a  large  impact  on  the  amount  of  dissolved  gases 
in  a  river  system.  Even  though  the  water  passes  over  the  structure  for  only  a 
short  time,  the  water  flowing  over  a  spUlway  or  weir  entrains  air  bubbles, 
creating  significantly  more  air-water  surface  area  for  gas  transfer.  In  addition, 
the  high  turbulence  that  occurs  at  most  hydraulic  structures  will  increase  the 
transfer  rate  coefficients.  The  same  quantity  of  gas  transfer  that  normally 
would  occur  in  several  miles  in  a  river  can  occur  at  a  hydraulic  structure. 

The  transfer  of  oxygen  from  the  atmosphere  to  the  water  is  often  of 
interest;  therefore,  it  seems  logical  to  directly  measure  oxygen  transfer.  How¬ 
ever,  there  are  some  problems  associated  with  the  measurement  of  dissolved 
oxygen  (DO)  concentration.  If  the  DO  level  is  close  to  saturation  (within 
approximately  2.5  mg/®,  the  tremendous  imcertainty  associated  with  the  cur¬ 
rent  measurement  techniques  makes  the  estimation  of  gas  transfer  useless 
(Gulliver  and  Wilhelms  1992).  Also,  if  the  reservoir  is  stratified,  it  is  difficult 
to  predict  withdrawal  from  the  various  layers  with  the  required  precision,  and 
usually  impossible  to  sample  at  the  spillway  crest  (Gulliver  and  Rindels  1993). 
Because  the  required  field  conditions  for  accurate  DO  measurement  often  do 
not  occur,  other  measurement  techniques,  such  as  the  tracer  technique,  are 
used. 

The  first  tracer  technique  for  reaeration  rate  measurement  in  open  channels 
was  introduced  by  Tsivo^ou  (1968).  The  basic  assumption  behind 
Tsivoglou’s  technique  is  that  the  ratio  between  the  desorption  of  the  tracer  gas 
and  the  absorption  rate  of  oxygen  is  a  constant.  The  ratio  of  tracer  gas  transfer 
rate  to  the  oxygen  transfer  rate  can  be  determined  through  laboratory  experi¬ 
ments  involving  the  simultaneous  transfer  of  both  gases.  Since  the  driving 
force  of  gas  transfer  is  the  partial  pressure  difference  between  the  concentra¬ 
tion  in  equilibrium  with  the  atmosphere,  which  is  linearly  proportional  to  the 
partial  pressure  in  the  atmosphere  and  that  dissolved  in  water,  tracers  that  have 
no  significant  atmospheric  component  should  be  chosen.  Thus,  measurements 
can  be  made  with  a  tracer  saturation  concentration  of  essentially  zero. 

Tsivoglou  injected  radioactive  kiypton-85  into  a  river  and  measured  its  desorp¬ 
tion  along  the  river  reach.  Tsivo^ou’s  technique  worked  rather  well  but 
required  the  use  of  radioactive  materials. 
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Rathbun  et  al.  (1978)  modified  Tsivoglou’s  method  by  substituting  propane 
and  ethylene  in  place  of  krypton-85  to  avoid  the  use  of  radioactive  materials. 
The  assumptions  made  were  similar  to  Tsivoglou’s  with  changes  made  for  the 
desorption  rates  of  propane  and  ethylene.  One  major  restriction  is  that  mea¬ 
surable  quantities  of  the  tracer  could  only  be  obtained  by  bubbling  the  gas 
through  pneumatic  diffusers,  requiring  additional  field  equipment.  A  more 
extensive  review  of  gas  tracer  techniques  is  given  in  McCutcheon  (1989). 

Thene  and  Gulliver  (1990)  developed  a  headspace  measurement  technique 
for  propane  gas  tracer  and  subsequently  found  measurable  amounts  of  naturally 
occurring  methane  in  the  samples.  The  unanticipated  quantities  of  methane 
held  promise  for  the  use  of  methane  as  a  tracer  gas.  Methane  is  produced  in 
the  sediments  as  a  by-product  of  the  anaerobic  decomposition  of  organic  mate¬ 
rials.  No  injection  of  artificial  tracer  gases  is  needed,  which  can  be  a  substan¬ 
tial  savings  in  cost  and  effort  at  the  larger  hydraulic  structures.  In  addition, 
the  uncertainty  associated  with  the  transverse  mixing  of  the  injected  tracers  is 
eliminated.  The  technique  is  more  effective  at  low-head  structures  since  strati¬ 
fication  of  methane  has  been  observed  at  higher  head  structures  (Wilhelms  et 
al.  1993). 

McDonald  and  Gulliver  (1992)  measured  oxygen  and  methane  transfer  over 
several  hydraulic  structures.  After  correction  for  diffusivities,  oxygen  and 
methane  transfer  efficiencies  were  foimd  to  be  comparable  at  a  given  structure 
except  when  the  entrained  air  bubbles  were  pulled  to  a  depth  in  the  tail  water. 
Because  of  the  high  static  and  dynamic  pressures  that  a  bubble  experiences  in 
a  plunging  jet  flow,  the  concentration  of  atmospheric  gases  will  be  greater  in 
the  bubble  than  at  the  water  surface.  The  rate  of  gas  transfer  will  correspond¬ 
ingly  be  higher  for  oxygen  and  nitrogen  due  to  the  higher  saturation  concen¬ 
tration  of  the  water  in  contact  with  the  bubble.  This  problem  was  delineated 
by  determining  an  “effective”  saturation  concentration  finm  oxygen  and 
methane  measurements  that  incorporates  the  higher  pressure  exerted  on  bubbles 
in  the  plunge  pool. 

This  report  describes  the  development  and  testing  of  the  methane  tracer 
technique  for  determining  gas  transfer  rates  at  low-head  hydraulic  structures. 
The  effect  of  the  plunge  pool  on  the  gas  transfer  rate  is  also  examined  by 
comparing  the  tracer  gas  transfer  rate  with  the  rate  for  oxygen  and  nitrogen. 
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2  Methane  in  the  Water 
Column 


Process  of  Methane  Production 

Methane  is  produced  as  a  by-product  of  the  anaerobic  decomposition  of 
organic  materid.  Methanogenesis  is  the  terminal  process  in  a  chain  of 
decomposition  processes  in  the  anaerobic  hypolimnion  and  represents  a  major 
mechanism  by  which  carbon  leaves  the  sediments  (Strayer  and  Tiedje  1978). 
Methane  transfer  out  of  the  sediments  and  hypolimnia  is  by  vertical  diffusion 
and,  if  enough  methane  is  produced,  by  bubble  ebullition. 

Methane  production  is  a  two-stage  process.  In  the  first  stage,  a  group  of 
faculative  and  obligate  bacteria,  termed  acid  formers,  turn  proteins,  carbohy¬ 
drates,  and  fats  into  fatty  acids  by  hydrolysis  and  fermentation.  Methane- 
producing  bacteria  then  convert  these  acids  into  methane.  Some  alcohols  from 
carbohydrate  fermentation  can  also  be  converted  to  methane  by  methane- 
producing  bacteria  (Wetzel  1975). 

The  concentration  of  acetate,  Hj,  and  CO2,  the  major  substrates  for  methane 
production,  is  highly  variable  and  is  dependent  upon  the  organic  material. 
Acetate  is  the  preferred  substrate  for  methanogenesis  at  low  partial  pressures  of 
H2.  As  the  partial  pressure  of  hydrogen  rises,  bicarbonate  is  the  preferred 
substrate. 

The  reaction  of  methane  production  from  the  reduction  of  carbon  dioxide  is 
given  as: 

CO2  +  8//^  ^  C//4  +  2Hp  (1) 

In  the  second  process,  the  production  of  methane  from  acetic  acid  can  be 

shown  as: 

CH^COOH  CH^  +  CO2  (2) 
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The  later  process  is  thought  to  produce  about  70  percent  of  methane  in  the 
water  column  with  much  of  the  remainder  coming  from  the  reduction  of  COj. 
Fermentation  of  other  acids  are  of  less  importance. 

Bacteria  that  produce  methane  are  strictly  anaerobic.  Even  slight  traces  of 
oxygen  are  toxic  to  these  bacteria.  They  consist  of  four  major  genre;  rod¬ 
shaped,  nonsporiferous  Methanobacterium,  rod-shaped  sporiferous 
MethanobaciUus,  the  spherical  Methanococcus,  and  Methanosurcina.  These 
bacteria  have  a  wide  temperature  range  for  growth.  Zeikus  and  Winfrey 
(1976)  found  methanogenesis  occurring  at  temperatures  from  4  to  45  °C, 
although  the  rates  were  highly  temperature  dependent.  A  change  of  12  °C  in 
the  temperature  of  the  sediments  during  the  change  of  seasons  was  associated 
with  a  1(X)-  to  400-fold  increase  in  the  rate  of  methanogenesis. 

Methanogenesis  is  severely  limited  by  sulfate  reduction.  Waters  that  have 
high  methane  concentration  generally  have  low  sulfate  concentrations  and  vice 
versa.  This  is  because  the  sulfate-reducing  bacteria  successfully  compete  with 
the  methane-producing  bacteria  for  and  acetate,  both  precursors  of  methane 
production  (Abram  and  Nedwell  1978).  The  reduction  of  sulfate  in  the 
sediments  is  described  by 

so;  +9H*  HS-  + 

MiUimolar  quantities  of  sulfate  will  inhibit  methane  production.  The  sulfate- 
reducing  bacteria  have  a  higher  affinity  for  hydrogen  and  acetate  than 
methane-producing  bacteria,  thus  keeping  the  pool  of  these  substrates  at  a  level 
too  low  to  be  used  by  the  methane-producing  bacteria  (Lovley  et  al.  1982).  In 
addition,  thermodynamic  calculations  can  be  used  to  predict  the  exclusion  of 
methane  production  in  sulfate  containing  sediments. 

Very  little  of  the  methane  produced  in  lakes  or  reservoirs  escapes  to  the 
atmosphere  because  of  the  relatively  quiescent  nature  of  the  water  body  or  the 
low  turbulence  intensities  near  the  water  surface.  Rudd  (1979)  measured  eva¬ 
sion  rates  of  5  percent  of  the  total  methane  produced.  Fallon  et  al.  (1980) 
measured  an  evasion  rate  of  9  percent  in  Lake  Mendota,  Wisconsin.  The  small 
evasion  rates  are  due  to  the  presence  of  methane-oxidizing  bacteria,  with  aero¬ 
bic  methane  oxidation  that  is  given  by: 

5CH^  +  Wj  2[CH;  +  3CO2  +  8//2O  (4) 


Because  the  residence  time  over  a  hydraulic  structure  is  short,  it  is  assumed 
that  no  significant  methane  oxidation  occurs  over  the  time  it  takes  a  methane 
molecule  to  pass  through  a  structure.  For  example,  while  the  half-life  of  meth¬ 
ane  over  a  structure  is  on  the  order  of  seconds  due  to  the  large  surface  area 
and  high  turbulence  of  the  flow,  the  half-life  of  methane  in  the  river  and  reser¬ 
voir  is  on  the  order  of  hours.  Therefore,  it  is  assumed  that  the  only  pathway 
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for  a  methane  molecule  to  leave  the  water  column  during  the  retention  time 
over  a  structure  is  by  evasion  to  bubbles  or  the  atmosphere. 


Methane  Concentrations 

Methane  was  sampled  in  a  number  of  reservoirs  to  estimate  how  often  one 
could  expect  relatively  high  (above  the  limit  of  quantification,  LOQ)  methane 
concentrations.  Methane  concentrations  were  sampled  in  seven  reservoirs 
during  fall,  winter,  and  spring;  in  all  but  one  reservoir,  they  were  weU  above 
the  current  LOQ  of  0.3  pg/«  (the  LOQ  is  developed  in  Chapter  6).  Sample 
methane  isopleths  obtained  from  the  Coon  Rapids  Dam  are  shown  in  Figures  1 
and  2.  The  concentrations  varied  from  119  to  4  pg/f  in  the  six  reservoirs.  AU 
six  reservoirs  were  fairly  shallow,  with  depths  ranging  from  3  m  (10  ft)  to  9  m 
(30  ft).  The  exception  was  one  relatively  deep  (18-m  or  60-ft)  reservoir  with  a 
methane  concentration  of  approximately  0.5  pg/C.  This  reservoir  was  held 
back  by  a  hollow  gravity  dam,  and  the  leakage  into  the  hollow  center  had  a 
strong  hydrogen  sulfide  smell.  As  mentioned  previously,  methanogenesis  is 
severely  limited  by  sulfate  production,  which  is  believed  to  be  occurring  in  this 
reservoir. 

Methane  was  supersaturated  at  aU  locations  sampled.  This  agrees  with 
results  from  de  Angelis  and  LiUey  (1987).  The  bacterial  oxidation  of  methane 
and  methane  evasion  rates  (i.e.,  the  pathways  by  which  methane  escapes  from 
the  water  column)  ate  inadequate  to  remove  all  the  methane  from  these 
reservoirs. 

Methane  concentrations  were  also  not  highly  stratified  in  the  shallow  reser¬ 
voirs  except  under  ice  cover,  where  methane  was  found  at  many  sites  at  much 
larger  concentrations  just  under  the  ice  than  near  the  sediments.  This  may  be 
due  to  methane  bubbles  rising  from  the  sediments.  The  ice  cover  would  block 
methane  evasion  to  the  atmosphere,  leaving  the  water  just  under  the  ice  at  a 
higher  concentration.  Another  thought  is  that  perhaps  the  ice  itself  contains  a 
high  amount  of  methane,  and  the  water  in  contact  with  the  melting  ice  in  early 
spring  will  also  have  a  correspondingly  high  methane  concentration.  Methane 
stratification  may  also  be  due  to  inputs  into  the  river/reservoir. 

Results  of  these  field  measurements  show  that  under  conditions  of  thermal 
stratification  in  reservoirs,  methane  will  also  be  stratified,  although  not  as 
strongly  as  oxygen.  The  highly  active  photic  zone  near  the  reservoir  surface 
will  enhance  the  naturally  occurring  oxygen  stratification.  Such  is  not  the  case 
for  methane.  The  primary  service  that  methane  can  provide  as  a  gas  tracer, 
however,  would  be  as  a  substitute  for  O2  and  N2  transfer  at  structures  associ¬ 
ated  with  nonstratified  reservoirs,  when  the  two  gases  are  close  to  saturation. 
Methane  will  likely  be  present  at  measurable  quantities  and  will  be  well  above 
saturation. 
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Figure  1.  Methane  isopleth  from  Coon  Rapids  Dam,  01  August  1989 
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Chapter  2  Methane  in  the  Water  Column 


Figure  2.  Methane  isopleth  from  Coon  Rapids  Dam,  06  October  1989 


3  Theory  and  Analysis 


Gas  Transfer  at  Hydraulic  Structures 


Gas  transfer  is  typically  described  as  a  first-order  process  in  which  the  rate 
of  change  of  the  gas  concentration  in  the  water  is  linearly  dependent  on  the 
concentration  difference  or  driving  force.  The  flux  of  any  volatile  dissolved 
chemical  across  an  air-water  interface  may  be  given  by  the  equation: 


VdC 
A  dt 


f  ^ 


(5) 


where 

F  =  mass  flux  rate  per  unit  surface  area 
V  =  Control  volume  over  which  C  and  A  are  measured 
A  =  Air-water  interface  surface  area 
C  =  Concentration  of  dissolved  chemical  in  water 
t  =  Time 

Ki  =  Liquid  film  coefficient 

C„  =  Concentration  of  chemical  in  air 

H  =  Henry’s  law  constant,  an  equilibrium  partitioning  constant 

C JH  is  often  called  the  saturation  concentration,  C^.  Because  the  chemicals 
oxygen,  nitrogen,  and  gas  tracers  of  interest  are  volatile,  only  the  liquid  phase 
resistance  to  transfer  (Kl)  is  considered.  The  gas  phase  resistance  is  much  less 
than  the  liquid  phase  resistance,  and  is  therefore  neglected. 
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For  transfer  at  a  hydraulic  structure,  Equation  5  can  be  integrated  from 
upstream  of  the  structure  to  downstream  of  the  structure,  resulting  in  an  equa¬ 
tion  for  transfer  efficiency  (Gulliver  and  Rindels  1993): 


E  =  £££ll=l-expj 


-  \Kj^t 


-J* 


(6) 


where  the  subscripts  “m”  and  “d”  refer  to  measurement  locations  upstream  and 
downstream  of  the  structure,  respectively.  A  transfer  efficiency  of  1.0  means 
the  full  gas  transfer  up  to  the  saturation  value  has  occurred.  No  gas  transfer 
would  correspond  to  E  =  0.0.  A  transfer  efficiency  of  greater  than  1.0  means 
the  gas  has  become  supersaturated  downstream.  For  gases  that  do  not  have 
appreciable  concentrations  in  the  atmosphere  Gike  methane),  Q  and  Q  are 
close  to  zero,  and  transfer  efficiency  can  be  written  as; 


E  = 


(nonatmospheric  gas  tranter) 


Indexing  Gas  Transfer 


(7) 


Gulliver  et  al.  (1990)  developed  an  indexing  relationship  to  directly  com¬ 
pare  liquid  phase  controlled  transfer  efficiencies,  such  as  those  considered 
herein,  at  hydraulic  structures.  This  relationship  can  be  used  to  index  transfer 
measurements  to  a  common  temperature  and  a  common  gas  (most  often  oxy¬ 
gen  is  used  in  gas  transfer  measurements).  The  relationship  is  given  as: 


where  E^  is  the  measured  transfer  efficiency,  i.e.,  the  transfer  efficiency  of 
methane,  and  is  the  equivalent  transfer  efficiency  of  the  index  gas  ,  i.e., 
oxygen,  at  the  index  temperature.  The  parameter/;  is  given  as 


fi  -fgfi 


(9) 


where /g  is  the  diffusivity  correction  and/j  is  the  temperature  correction, 
approximated  by  a  curve  fit  to  the  original  equation  as; 

/j.  =  1.0  +  0.02103(1-20)  +  8.261x10-5(1-20)2  (10) 

where  T  is  the  temperature  in  °C.  If  the  methane  transfer  measurements  are 
made  at  the  index  temperature,  no  temperature  indexing  is  needed  and/,  is  set 
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equal  to  1.  The  parameter/^  may  be  determined  from  the  parametric 
relationship  given  by  Hanratty  (1991)  and  Rathbun  (1990)  to  be: 


(11) 


where  D  is  the  diffusivity  of  the  measured  compound,  methane,  and  is  the 
diffusivity  of  oxygen.  Rathbun  (1990)  and  Gulliver  et  al.  (1990)  give  tech¬ 
niques  to  estimate  D  for  various  compounds.  The  diffusivity  of  methane  in 
water  was  researched  by  Thene  and  Gulliver  (1990),  who  searched  Chemical 
Abstracts  finom  1965  through  1986  and  found  14  measured  diffusivities  with  a 
mean  value  of  1.69  x  10"’  m^/sec  at  20  °C.  The  diffusivity  of  oxygen  (2.18  x 
10'^  m^/sec)  and  Equation  1 1  were  then  used  to  compute  an  value  of  0.88 
for  the  indexing  of  methane  transfer  measurements  to  oxygen  and  vice  versa. 


Comparison  of  Oxygen  and  Methane  Transfer 
Efficiencies 

Simultaneous  oxygen  and  methane  measurements  aUow  for  the  comparison 
of  their  transfer  efficiencies.  Figure  3  shows  the  measured  oxygen  transfer 
efficiencies  plotted  against  the  methane  transfer  efficiencies  indexed  to  oxygen 
for  five  low-head  structures.  The  oxygen  transfer  efficiency  measurement  was 
generally  equal  to  or  higher  than  the  methane  transfer  efficiency,  after  indexing 
to  an  equivalent  oxygen  transfer,  indicating  a  discrepancy  between  the  two 
simultaneous  measurements. 

One  possible  cause  of  this  discrepancy  could  be  that  the  indexing  between 
methane  and  oxygen  transfer  was  incorrect.  There  is  some  debate  over  the 
power  to  be  applied  to  the  diffusivity  ratio,  which  some  have  said  should  be 
two-thirds  rather  than  one-half  at  low  turbulence  levels.  Aside  from  the  fact 
that  flow  at  prototype  hydraulic  structures  is  highly  turbulent,  this  difference  in 
power  would  not  have  a  significant  effect  on /^.  The  use  of  an  incorrect^ 
could  not  account  for  the  variations  in  inequity  between  indexed  methane  and 
measured  oxygen  transfer  efficiencies  given  in  Figure  3.  It  would  simply  bias 
the  results  to  one  side  or  the  other  of  perfect  agreement. 

A  more  likely  cause  of  the  incongruency  of  the  methane  and  oxygen  trans¬ 
fer  efficiencies  is  due  to  the  effect  of  tail  water  on  gas  transfer.  As  illustrated 
in  Figure  4,  the  plunging  water  jet  at  the  surface  will  carry  entrained  air  bub¬ 
bles  to  some  depth.  The  hydrostatic  force  of  the  water  column  will  result  in  a 
greater  pressure  inside  these  entrained  bubbles  as  the  bubbles  are  pulled 
deeper.  Additionally,  as  the  jet  reaches  the  bottom  of  the  plunge  pool,  some 
force  acts  on  the  flow  to  remove  the  vertical  component  of  velocity,  deflecting 
the  spillway  jet  horizontally  out  of  the  plunge  pool.  This  force  also  increases 
the  pressure  on  the  entrained  bubbles.  The  increased  pressures  yield  a  greater 
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Figure  4.  Air  entrainment  at  plunge  pool  surface 
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partial  pressure  of  oxygen  in  the  bubbles,  and  hence,  an  increased  saturation 
concentration  than  would  occur  at  atmospheric  pressure. 

This  phenomenon  increases  the  transfer  of  atmospheric  gases  such  as  oxy¬ 
gen  and  nitrogen,  since  the  bubbles  will  contain  significant  amounts  of  these 
gases.  However,  for  methane  and  other  gases  that  do  not  have  sigmficant 
atmospheric  concentrations,  the  hydrostatic  and  dynamic  forces  acting  in  the 
plunge  pool  will  not  increase  the  gas  transfer  rate  since  the  partial  pressures  of 
these  gases  in  a  bubble  will  be  close  to  zero  in  most  cases.  This  will  be 
shown  for  one  typical  structure  later  in  the  report. 

The  integration  of  Equation  5  to  result  in  Equation  6  assumed  that  the 
saturation  concentration  was  constant.  For  methane  transfer,  the  saturation 
concentration  is  essentially  zero  regardless  of  the  pressure  within  the  bubble, 
and  therefore  Equation  6  is  applicable.  However,  for  atmospheric  gases,  this  is 
obviously  not  true.  Thus,  an  oxygen  transfer  efficiency  measurement  should 
assume  a  saturation  concentration  that  corresponds  to  pressures  somewhat 
higher  than  atmospheric.  Determination  of  this  pressure  and  corresponding 
saturation  concentration  is  one  of  two  primary  objectives  of  this  report  The 
primary  justification  is  that  when  a  tracer,  such  as  methane,  is  being  used  to 
predict  an  oxygen  or  nitrogen  transfer  efficiency,  or  when  the  oxygen  or  nitro¬ 
gen  transfer  efficiency  is  being  used  to  predict  the  transfer  of  other  non- 
atmospheric  gases,  substantial  errors  could  result.  In  addition,  accurate 
prediction  of  oxygen  and  nitrogen  transfer  at  hydraulic  structures  with  a  plunge 
pool  is  not  possible  without  some  quantification  of  this  tail  water  effect 
(Johnson  1984).  A  technique  is  therefore  needed  to  adjust  for  the  higher  pres¬ 
sure  that  the  entrained  bubbles  experience  in  the  tail  water  pool  of  a  hydraulic 
structure. 


Effective  Saturation  Concentration  and  Effective 
Bubbie  Depth 

The  technique  used  to  adjust  for  the  higher  pressure  that  the  entrained  bub¬ 
bles  experience  in  the  plunge  pool  assumes  that  methane  and  either  oxygen  or 
nitrogen  concentration  measurements  have  been  made  at  the  structure.  The 
process  is  then  to  (a)  measure  the  methane  transfer  efficiency,  E„,  of  the  struc¬ 
ture  using  Equation  7;  (b)  then  index  the  methane  transfer  efficiency  to  an 
oxygen  or  nitrogen  transfer  efficiency  with  Equation  8;  and  (c)  then  use  the 
upstream  and  downstream  dissolved  oxygen  (DO)  or  dissolved  nitrogen  mea¬ 
surements  with  the  indexed  transfer  efficiency  to  compute  an  “effective”  satu¬ 
ration  concentration  for  oxygen  or  nitrogen, 


C 


St 


C,-Cu 

EiO 


*  c. 


(12) 
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where  Eio  is  the  measured  methane  transfer  efficiency  indexed  to  an  oxygen  or 
nitrogen’ transfer  efficiency  and  C„  and  Q  are  the  measured  DO  or  dissolved 
nitrogen  concentrations  upstream  and  downstream,  respectively. 

The  effective  saturation  concentration  is  therefore  a  mean  of  the  oxygen  or 
nitrogen  saturation  concentrations  that  the  bubbles  experience,  weighted  by 
their  impact  on  gas  transfer.  This  weighted  mean  saturation  concentration  is 
not  determined  by  following  bubble  paths  through  a  plunge  pool,  but  is  simply 
inferred  from  simultaneous  methane  and  oxygen  or  nitrogen  concentration 
measurements. 

The  effective  saturation  concentration  can  be  expressed  as  an  “effective 
depth,”  or  the  depth  at  which  a  bubble  would  be  held  under  hydrostatic  pres¬ 
sure  to  have  the  same  equilibrium  concentration  as  the  effective  saturation 
concentration.  This  depth  could  be  compared  with  tail  water  depth,  specific 
discharge,  head,  etc.,  at  given  structures.  It  will  be  used  for  comparison  with 
these  parameters  because  of  the  perspective  it  gives  relative  to  tail  water  depth. 
The  effective  depth  is  computed  as: 


r 


-  ^p\ 


_p 


comp 


am 


(13) 


where 

Dgff  =  Effective  bubble  depth,  m  or  ft 

K  =  Conversion  constant  from  atmosphere  to  depth  of  water  (=33.9  ft 
HjO/atm  or  10.34  m  HjO/atm),  assuming  hydrostatic  pressure 

=  Effective  saturation  concentration  of  oxygen  or  nitrogen,  mg/4, 
from  Equation  12 

H  =  Henry’s  Law  coefficient  for  oxygen  or  nitrogen  (atm  4/mg) 

comp  =  Fraction  of  gas  in  atmosphere  (0.2095  for  oxygen  and  0.7808  for 
nitrogen) 

^atm  -  Barometric  pressure,  atm 

Impact  of  Assuming  a  Zero  Saturation 
Concentration  for  Methane 

Equation  7  assumes  the  saturation  concentration  of  methane  is  small  since 
the  concentration  of  methane  in  the  atmosphere  is  essentiaUy  zero.  However, 
the  concentration  of  methane  in  a  bubble  may  not  approximate  zero  due  to 
methane  transferring  into  the  bubble  from  the  water  as  the  bubble  is  dragged 
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through  the  plunge  pool.  Assuming  all  methane  transfer  occurs  across  a  bub¬ 
ble  interface,  a  mass  balance  of  methane  across  the  hydraulic  structure  is  writ¬ 
ten  as: 


(14) 


where 

C  =  Concentration  of  methane  in  a  bubble  as  it  is  released  to  atmosphere 
=  Unit  discharge  of  water,  m^/sec 
=  Unit  entrainment  of  air,  m^/sec 

C„  =  Concentration  of  methane  in  water  upstream  of  stmcture 

Cj  =  Concentration  of  methane  in  water  downstream  of  structure 

Measurements  taken  at  the  Rum  River  gated  OG  spillway  were  used  to 
assess  the  impact  of  assuming  a  methane  saturation  concentration  of  zero. 

C„  and  Cj  were  measured  at  16.58  and  8.72  mg/m^(ppb),  respectively,  q^  was 
calculated  from  the  gate  opening,  head,  and  geometry  of  the  stmcture  to  be 
0.167  mVsec.  The  specific  discharge  of  the  entrained  air  was  measured 
directly  using  a  1-  by  1-m  square  hood  designed  and  built  by  Mr.  Steven  C. 
Wilhelms.  The  hood  was  placed  on  the  water  surface  in  the  plunge  pool  cap¬ 
turing  the  air  released  to  the  atmosphere  by  the  bubbles.  The  hood  constriaed 
into  a  2-cm-diam  pipe,  where  a  hot  film  anemometer  was  placed  to  measure 
velocity,  enabling  the  air  discharge  coming  out  of  the  water  to  be  calculated. 
These  measurements  are  given  in  the  Appendix  A.  The  hood  was  moved 
laterally  and  longitudinally  along  the  tailrace  in  the  region  where  air  bubbles 
were  escaping  and  the  air  flow  rates  were  summed  to  yield  the  total  air  dis¬ 
charge.  The  specific  air  discharge  was  then  calculated.  For  this  experiment, 
q^j  was  computed  from  these  measurements  to  be  0.125  m^/sec. 

Then,  using  Equation  14  to  compute  the  released  concentration  of  methane 
for  a  bubble: 


c  .  H£!(;n!^(I6.58-8.72)m*/m’-10.5ms/m’  (15) 

"  0.125(m2/sec) 


The  equilibrium  concentration  of  methane,  C^,  in  a  bubble  downstream  of  the 
Stmcture  can  be  found  from  the  Henry’s  Law  relationship: 


(16) 
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where  Hem  =  dimensionless  Henry’s  Law  constant  for  methane  =  17.98  at 
0  °C.  Using  Equation  16  to  calculate  C^. 


=  8.72  mg/m^  (17.98)  =  156.8  mg/m^ 


(17) 


Therefore,  the  concentration  of  methane  in  a  bubble  as  it  is  released  to  the 
atmosphere  is  6.7  percent  of  the  equilibrium  concentration  of  methane  in  that 
bubble. 

The  transfer  efficiency  of  methane  can  be  calculated  by  rewriting  Equa¬ 
tion  6  in  terms  of  Henry’s  Law  constant: 


(18) 


where  Cf,  is  an  average  value  of  the  methane  concentration  in  the  bubble,  or: 


C,  = 


1 


-  <7m 

\  J 


1  + 


D 


10.3 


(19) 


The  concentration  of  methane  in  the  bubble  varies  as  the  bubble  is  washed 
through  the  plunge  pool.  At  the  inception  of  the  bubble,  the  concentration  is 
approximately  zero,  since  there  is  relatively  little  methane  in  the  atmosphere. 
As  the  bubble  is  released  to  the  atmosphere  at  the  end  of  the  plunge  pool,  the 
concentration  is  C„.  The  first  term  in  Equation  19  is  an  arithmetic  mean  of 
the  end  point  concentrations.  The  second  tenn  of  Equation  19  accounts  for  the 
increase  in  concentration  due  to  the  reduction  in  volume  of  a  discrete  bubble 
as  it  is  dragged  to  some  depth  in  the  plunge  pool.  For  the  example  used 
herein,  was  computed  from  a  comparison  of  oxygen  and  methane  transfer 
measurements  to  be  0.8  m. 


Solving  Equations  18  and  19: 


= 


8.72-16.58 


(  5.7  ^ 

-  16.58 

[l7.98j 

=  0.48 


C,  =  1  (10.5) 


1  + 


0.832 

10.3 


=  5.7  mg/m  • 


(20) 

(21) 


When  the  concentration  of  methane  in  a  bubble  is  ignored,  Q  =  0  and 
Equation  7  becomes: 
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(22) 


8.72-16.58 
-  16.58 


0.47 


Thus,  the  difference  in  transfer  efficiencies  calculated  by  accounting  for  a 
concentration  of  methane  in  the  bubble  and  by  ignoring  the  concentration  of 
methane  in  the  bubble  is  only  1.9  percent  of  the  transfer  efficiency.  This  is 
well  within  the  uncertainty  of  the  measurements  and  was  incorporated  as  a  bias 
imcertainty  in  the  calculations. 

Further  calculations  using  a  predictive  relationship  for  the  (Ervine  and 
Elsawy  1975)  calibrated  to  these  measurements  for  entrained  air  discharge 
indicated  that  the  bias  of  roughly  2  percent  would  apply  to  all  discharges  at 
this  structure. 
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4  Sampling  Technique 


Oxygen  and  methane  measurements  were  taken  at  all  of  the  stnictures 
tested.  Additionally,  total  dissolved  gas  measurements  were  taken  at  the  Rum 
River  stracture  to  determine  dissolved  nitrogen  concentration.  At  each  site,  an 
upstream  sampling  location  was  chosen  as  close  to  the  structure  as  possible  to 
ensure  that  the  water  sampled  was  representative  of  the  water  passing  over  the 
stracture.  The  downstream  sampling  point  was  chosen  close  to  the  stracture, 
yet  downstream  of  the  air  entrained  or  bubbly  region  of  the  plunge  pool.  At 
both  upstream  and  downstream  sampling  points,  measurements  were  taken  near 
the  water  surface,  at  mid-depth,  and  near  the  bottom  to  check  for  stratification 
of  methane  or  oxygen. 


Dissolved  Oxygen  Sampling 


Dissolved  oxygen  measurements  were  performed  in  situ  with  a  YSI  No.  58 
dissolved  oxygen  probe  and  meter.  The  probe  was  lowered  to  the  appropriate 
sampling  depth  along  a  weighted  tether  to  the  probe  from  drifting  to  the 
surface  in  locations  of  high  stream  velocity.  The  meter  reading  of  either  con¬ 
centration  (mg/«)  or  percent  of  saturation  was  recorded.  The  saturation  con¬ 
centration,  Cj,  was  detennined  as; 

^  _  0.0295  (23) 

*  H 


where  H  is  Henry’s  Law  and  is  atmospheric  pressure.  Gulliver  and 
Rindels  (1986)  found  the  actual  saturation  concentration  in  river  water  to  be 
98  ±  2  percent  of  the  values  listed  in  the  literature. 


Total  Dissolved  Gas  Sampling 


The  total  dissolved  gas  concentration  was  measured  in  situ  at  the  Rum 
River  structure  using  a  modified  “gasometer”  (Bouck  1982)  consisting  of  a  gas 
peimeable  tubing  connected  to  a  pressure  transducer  as  shown  in  Photo  1.  The 
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gasometer  was  also  lowered  to  the  appropriate  depth  along  the  weighted  tether. 
After  allowing  the  gasometer  to  equilibrate,  the  voltmeter  reading  from  the 
pressure  transducer  was  recorded.  In  most  cases,  it  took  approximately  15  to 
20  min  for  the  gasometer  to  equilibrate.  The  transducer  was  calibrated  before 
each  field  trip  to  yield  a  known  pressure  for  a  given  signal,  as  shown  in  Fig¬ 
ure  5.  The  scatter  of  the  calibration  points  shown  in  Figure  5  is  within  the 
range  of  uncertainty  reported  by  the  manufacturer.  The  slope  of  the  calibration 
was  determined  from  linear  regression  and  the  uncertainties  reported  by  the 
manufacturer  were  included  in  the  uncertainty  of  our  measurements. 

The  dissolved  nitrogen  concentration,  Cf^^,  was  computed  as: 


where  is  the  Henry’s  Law  constant  for  nitrogen  and  partial 

pressure  of  nitrogen  inside  the  tubing  calculated  as: 

=  98.77%  (P  -  Fo,  -  \o)  (25) 

where  F  is  the  total  dissolved  gas  pressure  (measured  with  the  gasometer), 

Ph2o  vapor  pressure  of  water,  and  Pq^  is  the  partial  pressure  of  oxygen 
in  the  tubing  calculated  as: 

O^molc  fraction  in  dry  air  (26) 

where  O2  percent  is  the  percent  saturation  of  oxygen  in  the  water,  P^„  is  the 
barometric  pressure,  and  F;^  is  the  vapor  pressure  of  water.  The  mole  frac¬ 
tion  of  oxygen  in  dry  air  is  D.2095  (Weast  1976). 


Methane  Sampling 

Methane  sampling  was  performed  by  collecting  water  samples  in  the  field 
and  later  analyzing  those  samples  in  the  laboratory  using  the  headspace  analy¬ 
sis  technique  (described  in  Section  VI).  Water  samples  were  collected  in 
40-mL  borosilicate  glass  vials  with  Teflon-faced  septa  and  open-top  screw 
cj^s.  The  uncapped  empty  vials  were  loaded  in  the  sampling  device  shown  in 
Photo  2.  The  sampler  was  rapidly  lowered  to  the  appropriate  depth.  Water 
enters  the  sampler  through  the  side  tubings,  and  air  exits  the  sampler  through 
the  top  tube.  This  design  virtually  eliminates  bubbling  within  the  sampler  that 
could  result  in  stripping  of  methane  from  the  water.  Once  filled,  the  sampler 
was  removed  from  the  water.  The  vials  were  capped.  Teflon  side  down,  under 
water  while  stfll  in  the  sampler  in  order  to  reduce  the  volatilization  of  methane 
to  the  atmosphere.  The  capped  samples  were  packed  top  down  to  prevent  the 
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Figure  5.  Total  gas  meter  calibration 


loss  of  methane  due  to  bubble  formation  caused  by  temperature  change  and 
subsequent  leakage  through  the  septa  or  cap. 


Sample  preservation 

As  samples  were  not  analyzed  immediately  after  they  were  gathered,  some 
method  of  preservation  had  to  be  employed.  Schultz  et  al.  (1971)  concluded 
that  1  mL  of  37-percent  fonnalin  solution  in  a  60-mL  solution  preserved  the 
sample  for  7  days.  To  stop  bacterial  action,  other  researchers  have  raised  the 
pH  to  11,  added  HgQ  or  stored  samples  on  ice.  It  was  decided  to  test  0.5  and 
0.25  mL  of  37-percent  formalin  solution  to  preserve  the  samples.  Samples 
were  gathered  from  the  Anoka  Dam  on  the  Rum  River  and  from  the  Missis¬ 
sippi  River  in  Minneapolis.  As  can  be  seen  from  Figure  6,  the  methane  con¬ 
centrations  with  1/4-mL  formalin  and  with  no  formalin  decreased  with  time. 
The  samples  with  1/2-mL  formalin  were  relatively  constant.  Therefore,  it  was 
decided  to  use  1/2  mL  of  formalin  as  a  preservative  for  this  smdy.  The  forma¬ 
lin  was  injected  in  the  samples  as  soon  as  possible  after  collection  and  before 
transport. 


Figure  6.  Effect  of  formalin  on  methane  concentration 
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Sample  transporting 

The  samples  were  transported  to  the  gas  chromatograph  (GC)  laboratory  at 
the  University  of  Minnesota,  Department  of  Civil  and  Mineral  Engineering  for 
analysis.  Samples  obtained  near  the  University  were  transported  immediately 
by  car,  whereas  samples  obtained  elsewhere  were  shipped  via  overnight 
express  couriers. 

All  water  samples  collected  at  a  given  location  exhibited  a  certain  degree  of 
precision  uncertainty  in  methane  concentration.  However,  those  samples 
shipped  by  overnight  couriers  exhibited  such  a  large  precision  uncertainty  in 
methane  concentration  that  for  all  practical  purposes  the  samples  were  useless. 
Tables  1  and  2  show  that  sample  vials  collected  at  the  same  location  and  time, 
and  even  from  the  side  by  side  bottles  in  the  same  sampler  grab,  varied  in 
methane  concentrations  by  as  much  as  69  percent  While  there  is  a  possibility 
that  this  is  truly  representative  of  the  actui  concentrations,  it  is  extremely 
unlikely  that  there  is  that  much  of  a  variation  in  methane  concentrations  due  to 
the  turbulent  mixing  at  structures.  A  more  likely  explanation  for  the  large 
variances  in  concentrations  is  that  there  was  a  depletion  of  methane  from  some 
of  the  sample  vials  during  shipping. 

Although  the  samples  were  packed  upside  down,  there  is  no  way  of  know¬ 
ing  if  the  samples  remained  upside  down  for  their  entire  journey  to  Minneapo¬ 
lis.  Perhaps  bubbles  of  high  methane  concentration  formed  and  leaked  out  the 
bottle-septa  seal  or  through  a  piercing  hole  in  the  septa  while  inverted.  The 
potential  for  methane  loss  during  transport  is  increased  by  the  atmospheric 
pressure  drop  exerted  on  the  samples  during  airplane  transport.  Figure  7 
shows  possible  methane  evasion  routes.  In  light  of  the  methane  concentration 
randomness  of  the  samples  transported  by  air  freight,  an  alternative  method  of 
transport  is  clearly  needed. 


Alternative  sample  shipping  technique 

An  alternative  packing  and  shipping  technique  has  been  developed  and  is 
currently  being  tested.  The  technique  is  referred  to  as  the  “double  redundancy 
technique.”  The  basis  of  the  technique  is  that  the  sample  vials  are  packed 
upside  down  immersed  in  a  water  bath.  The  technique  is  doubly  redundant 
since  first  the  risk  of  the  of  the  inverted  sample  being  righted  is  reduced.  The 
package  handlers  will  be  less  likely  to  tip  over  a  package  labeled  “T  LIQUID- 
THIS  END  UP  T”  in  order  to  avoid  a  spill  (actually,  a  spill  would  be  very 
unlikely  due  to  the  watertight  packaging).  Secondly,  even  if  the  sample  vials 
are  righted,  there  is  less  potential  for  methane  escape  due  to  the  vials  being 
packed  in  a  water  bath.  Figure  8  shows  the  restricted  flow  through  a  potential 
leak  in  a  septa.  The  surface  tension  of  the  water  on  the  outside  face  of  a  septa 
will  prohibit  the  flow  of  gas  from  an  air  bubble  on  the  bottle  side  of  a  septa. 
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Table  1  (Concluded) 

Location 

Gate 

Opening 

ft 

Sample 

Depth 

ft 

Time  of 
Sampling 

Bottle 

No. 

[CHJ 

ng/« 

Downstream 

1.5 

5 

13:55 

62 

4.57  ±0.19 

2.0 

5 

14:10 

60 

4,67  ±0.13 

68 

5.30  ±  0.48 

63 

5.37  ±0.15 

61 

5.51  ±  1.34 

Table  2 

Methane  Concentrations 
September  23, 1992 

i,  Smithland  Lock  and  Dam, 

Location 

Gate 

Opening 

ft 

Sample 

Depth 

ft 

Time  of 
Sampling 

Bottle 

No. 

[CHJ’ 

Upstream 

2.0 

15 

11:57 

30 

3.08  ±  0.09 

22 

2.89  ±0.19 

41 

3.58  ±  0.29 

49 

2.75  +  0.14 

Downstream 

0.5 

15 

9:40 

9 

10.22  +  0.54 

10 

7.66  ±  0.55 

12 

12.16  ±0.49 

1.0 

15 

10:05 

6 

7.17  ±0.08 

14 

6.63  ±  0.75 

7 

6.95  ±  0.20 

16 

7.33  ±1.32 

1.5 

15 

10:25 

13 

5.55  ±  0.59 

25 

6.06  ±  0.87 

5 

6.71  ±  0.47 

2.0 

15 

10:55 

20 

6.74  ±  0.56 

21 

6.58  ±  0.56 

11:00 

35 

5.14  ±0.28 

28 

5.44  ±  0.30 

^  Uncertainties  include  only  precision  uncertainty. 
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possible  leakage 
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Potential  evasion  routes  of  methane  during  sample  transport 
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Figure  8.  Restricted  tlow  of  gas  from  sample  immersed  in  water  bath 


In  practice,  an  empty  1-gal^  paint  can  could  be  used  as  the  container.  The 
vials  could  be  cushioned  with  sponge  or  foam  packing.  The  packing  would 
also  serve  as  a  spacer  to  keep  the  septa  of  all  vials  below  the  water  surface 
inside  the  can,  regardless  of  the  upward  orientation  of  the  can.  No  results  on 
the  effectiveness  of  the  technique  are  available  as  of  this  writing  since  the 
technique  is  currently  being  tested. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units  is 
presented  on  page  viii. 
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5  Headspace  Analysis 
Technique 


A  slightly  modified  version  of  the  headspace  analysis  technique  developed 
by  Thene  and  Gulliver  (1990)  was  used  to  determine  the  concentration  of 
methane  in  the  water  samples  by  measuring  the  methane  concentration  of  the 
air  in  equilibrium  with  the  water. 

Calculation  of  Methane  Concentration  in  Water 

The  concentration  in  the  headspace,  Q,  is  given  by  Thene  and  Gulliver 
(1990)  as; 


Wi 


(27) 


where 

C  =  Concentration  in  the  water  prior  to  the  creation  of  the  headspace 
R  =  Universal  gas  constant 
T  =  Water  temperature 
Af  =  Molar  mass  of  methane 
H  =  Henry’s  Law  constant  for  methane 
=  Volume  of  headspace  in  sample  bottle 
=  Volume  of  water  in  sample  bottle 
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The  exact  volumes  of  headspace  and  water  in  the  samples  used  in  Equation  27 
were  determined  by  weighing  each  bottle  before  creating  the  headspace,  after 
creating  the  headspace,  and  when  empty. 


Creating  the  Headspace 


A  headspace  was  created  in  each  sample  vial  at  the  GC  laboratory.  Photo  3 
shows  the  method  used  to  make  the  headspaces.  A  drainage  needle  open  to 
the  atmosphere  was  inserted  in  the  septa  of  the  upside  down  vial.  Approxi¬ 
mately  10  mL  of  nitrogen  gas  were  injected  through  the  septa  at  the  top  of  the 
vial,  forcing  water  out  through  the  drainage  needle.  Samples  were  then  shaken 
vigorously  for  60  sec  to  strip  the  methane  from  the  water  into  the  headspace. 


Gas  Chromatography 

The  gas  chromatograph  (GC)  used  was  a  Hewlett  Packard  5890A  equipped 
with  a  flame  ionization  detector  (FID),  a  strip  chart  recorder,  and  an  electronic 
integrator  and  a  4-ft-long  5A  60/80  molecular  sieve. 

Five  separate  200-pC  volumes  were  withdrawn  hum  the  headspace  of  each 
bottle  with  a  250-p{  gas-tight  syringe  and  then  injected  into  the  GC.  At  least 
three  injections  were  needed  to  accurately  compute  the  precision  uncertainty  of 
the  measurement:  five  was  chosen  as  a  matter  of  convenience.  The  pressure  in 
the  headspace  was  kept  nearly  constant  by  counter-injecting  200  }iC  of  water 
into  the  sample  prior  to  withdrawal  of  the  GC  injection. 

The  water  used  to  replace  the  sample  volume  must  not  contain  methane. 
Headspace  analysis  of  de-ionized  water,  tap  water,  and  Mega-Pure  water 
showed  that  only  the  Mega-Pure  water  sample  did  not  give  a  GC  response. 
Therefore,  only  Mega-Pure  water  was  used  to  replace  the  extracted  200-p8 
volumes.  The  injection  of  1  mL  of  Mega-Pure  water  does  slightly  affect  the 
methane  concentration  of  the  headspace.  However,  the  effect  is  small  due  to 
methane’s  high  volatility  (McDonald  and  Gulliver  1992;  Thene  and  Gulliver 
1990). 

Calibration  of  the  GC  was  perfoimed  by  injecting  different  volumes  of 
100  ppm  (by  volume)  methane  standard  and  plotting  the  results.  A  curve  was 
then  fit  through  the  calibration  points  to  obtain  an  equation  relating  area  count 
to  mass  of  methane,  as  shown  in  Figure  9.  This  mass  was  then  related  to 
concentration  in  the  original  sample  as  described  in  McDonald  and  Gulliver 
(1992).  Uncertainties  in  the  concentration  measurements  were  also  computed 
to  the  95-percent  confidence  interval.  The  confidence  interval  was  typically 
±1  to  2  percent  of  the  measured  concentration. 
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Figure  9.  GC  Calibration 


A  quality  control/quality  assurance  program  was  also  carried  out  to  ensure 
that  the  concentrations  in  the  water  column  were,  within  the  confidence  limits, 
correctly  determined  from  the  GC  readings.  Tests  were  performed  on  the 
length  of  time  the  samples  could  be  preserved,  the  quantity  of  formalin 
required  for  preservation,  the  techmque  of  sample  injection  into  the  GC,  the 
impact  of  headspace  pressure  change  due  to  sampling,  the  GC  response  of 
formalin,  the  reuse  of  septa,  and  the  effect  of  bubble  formation  in  the  vials 
(McDonald  and  Gulliver  1992). 


Injection  technique 

It  was  discovered  that  the  speed  in  which  the  sample  was  injected  into  the 
GC  led  to  different  responses.  To  test  the  effects  of  injection  speed,  different 
volumes  of  standard  were  injected  into  the  GC  at  different  speeds.  Three 
speeds  were  tested.  The  first  speed  was  to  inject  as  fast  as  possible.  The  total 
time  to  pierce  the  septa,  depress  the  plunger,  and  withdraw  the  syringe  was 
under  0.5  sec.  The  second  speed  was  a  3*sec  process.  One  second  was  used 
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to  pierce  the  septa  and  to  get  the  syringe  in  place.  Two  seconds  were  used  to 
depress  the  plunger,  then  the  syringe  was  withdrawn  quickly.  The  third  speed 
tested  was  similar  to  the  second;  only  4  sec  total  were  required  for  the  entire 
injection  process.  As  shown  in  Figure  10,  the  fast  injection  yields  a  different 
response  from  the  other  injection  speeds.  Thus,  if  the  injection  speed  is 
0.6  sec  instead  of  0.5  sec,  the  results  could  be  different.  The  other  two  plots 
are  virtually  identical,  indicating  that  the  results  were  very  nearly  insensitive  to 
injection  speed  for  the  3-  and  4-sec  injection.  Therefore,  it  was  decided  to  use 
the  3-sec  injection  for  samples.  It  is  thought  with  the  fast  injection  that  some 
of  the  sample  is  lost  through  the  pierced  septa  in  the  injection  port  of  the  GC 
or  through  the  plunger  seal  in  the  gas-tight  syringe. 


Figure  10.  Impact  of  injection  speed  on  GC  response 


Limit  of  detection/iimit  of  quantification 

The  limit  of  detection,  LOD,  is  defined  as  the  lowest  concentration  that  is 
statistically  different  from  the  blank.  The  limit  of  quantification,  LOQ,  is  the 
concentration  level  above  which  quantitative  results  can  be  obtained  with  con¬ 
fidence.  In  most  cases,  the  LOD  and  LOQ  are  defined  as  three  and  ten  times 
the  standard  deviation  of  the  response  to  blank  runs,  respectively  (American 
Chemical  Society  1983).  Experience  has  shown  that  there  is  no  methane 
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response  to  blank  runs.  Consequently,  a  different  procedure  was  required.  It 
was  decided  that  an  appropriate  LOD  and  LOQ  could  be  determined  by 
injecting  various  masses  of  the  standard  gas  each  10  times  into  the  GC  and 
calculating  the  standard  deviation  of  each  subset  of  injections.  Then  the  LOQ 
was  defined  as  the  point  where  ten  times  the  standard  deviation  of  a  set  of 
injections  equaled  the  average  area  count  of  the  set  of  injections.  Or  LOQ  =  C 
when  C  =  lOo,  where  C  is  the  mean  concentration  and  o  the  standard  devia¬ 
tion  in  concentration  measurements.  Similarly,  the  LOD  was  defined  as  the 
point  where  three  times  the  standard  deviation  equaled  the  average  area  count. 

Using  this  definition  and  the  information  gathered  from  the  injection  of 
standard  gas.  Table  3  was  prepared.  This  table  shows  the  LOQ  to  lie  between 
the  mass  contained  in  3  and  5  p{  of  standard  gas.  By  linear  interpolation 
between  these  two  points,  the  mass  of  LOQ  was  equal  to  0.000211  pg. 


Table  3 

Data  Used  for  LOQ  Analysis 

Volume,  pi 

Mass,  pg 

Mean  AC 

10  0 

lOo/mean 

Equivalent  Water 
Concentration, 

150 

0.009963 

35,045.1 

5,539.3 

0.16 

14.7 

100 

0.006642 

22,996.1 

7,456.6 

0.32 

9.5 

50 

0.003321 

12,192.9 

4,356.6 

0.36 

4.9 

10 

0.000664 

2,886.8 

1,658.2 

0.57 

0.98 

5 

0.000332 

1,891.4 

1,118.3 

0.59 

0.49 

3 

0.000199 

1,007,1 

1,047.5  ■ 

1.04 

0.29 

Assuming  a  250-pf  sample  injection,  the  headspace  concentration  would  be: 


C  =  =  0.84pg/« 

“  volume  250p  t 


(28) 


Since  the  volume  of  headspace  is  roughly  10  mf  and  the  volume  of  water  is 
roughly  30  m{,  the  pre-headspace  water  concentration  is  computed  from  Equa¬ 
tion  12  as: 


C  =  (0.04-K).33)C^  =  0.37C^  =  0.3pg/« 


(29) 
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This  is  the  value  for  LOQ.  LOD  is  defined  as  3/10  of  the  LOQ  such  that. 


LOD  =  ^  LOQ  =  0.1jig/« 


(30) 


Bubble  Formation 


Bubbles  may  be  formed  in  the  samples  as  the  samples  are  transported  from 
cold  temperatures  to  warmer  temperatures.  A  calculation  was  thus  made  to 
determine  if  the  bubbles  formed  would  significantly  affect  the  methane  con¬ 
centration  in  the  sample.  For  this  calculation,  a  3-mm-diam  bubble  and  a 
water  methane  concentration  of  C  =  20  pg/{  were  assumed.  Equation  27  then 
becomes: 


=  20pg/«  (0.03740.00036)  =  535  pg/«  (31) 


The  percentage  of  the  total  methane  in  the  sample  bottle  that  is  in  the  bubble, 
percent  CH4b„bbie.  Is: 


percent  CH^ 


CV 

100..— 

CV^ 


=  0.96  percent 


(32) 


The  bubble  would  remain  in  the  vial  upon  headspace  injection  and  would 
become  a  part  of  the  headspace.  If  10  mL  of  water  is  removed,  while  30  mL 
remain,  the  headspace  concentration  for  this  example  would  be  high  by 
1/4(0.96  percent)  =  0.24  percent  This  percent  difference  is  insignificant  com¬ 
pared  with  the  precision  of  the  GC  analytic  technique  for  methane. 
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6  Results 


The  results  of  simultaneous  methane-oxygen  transfer  at  four  stmctures  are 
summarized  in  Table  4.  Descriptions  of  the  structures  are  included  in  Appen¬ 
dix  B.  Figure  11  shows  the  effective  depth  plotted  versus  specific  discharge  at 
Rum  River  for  six  separate  field  trips.  Figure  12  is  a  similar  plot  for  the 
Sl  Qoud  structure  on  15  March  1990  and  05  April  1993.  The  calculated 
effective  depth  at  the  Anoka  structure  is  reproducible  as  a  function  of  specific 
flow  rates  for  a  variety  of  water  temperatures,  upstream  oxygen  deficits,  and 
heads.  The  effective  depths  at  the  St.  Qoud  differ  from  those  at  the  Rum 
River  stracture  for  comparable  specific  discharges,  perhaps  due  to  differences 
in  head,  tail  water  depth,  and  plunge  pool  design. 

Generally,  as  the  specific  discharge  increases,  the  effective  depth  increases. 
At  low  discharges,  the  momentum  of  the  plunging  jet  is  insufficient  to  carry 
the  bubbles  deep  into  the  pool,  and  the  eff^ective  saturation  concentration  is 
close  to  local  atmospheric  concentration.  At  high  discharges,  the  momentum 
of  the  plunging  jet  carries  the  entrained  bubbles  deeper  into  the  pool,  and  the 
effective  saturation  concentration  is  significantly  higher  than  that  computed 
from  local  atmospheric  concentration  because  of  the  higher  pressures  that  the 
bubbles  experience.  Negative  values  of  effective  depth  are  listed  in  Table  4 
for  the  lower  discharges  at  the  Kost  and  St.  Qoud  stmctures.  There  is  no 
physical  explanation  for  this.  The  actual  effective  depth  values  are  most  likely 
positive  and  the  measurements  are  in  error.  In  each  case  where  a  negative 
effective  depth  is  calculated,  there  is  a  greater  uncertainty  in  the  reported  effec¬ 
tive  depth.  Thus,  these  negative  values  are  small  relative  to  the  measurement 
uncertainty. 

As  a  verification  of  the  principles  presented  herein,  effective  depths  were 
calculated  from  dissolved  nitrogen  measurements  and  from  dissolved  oxygen 
for  a  range  of  discharges,  using  methane  as  the  tracer  gas  for  both.  Nitrogen 
and  oxygen  are  both  present  in  quantity  in  the  atmosphere  and  therefore  in  the 
entrained  bubbles.  As  a  bubble  is  puUed  to  some  depth  in  the  plunge  pool,  the 
same  hydrostatic  forces  increase  the  partial  pressures  of  both  oxygen  and  nitro¬ 
gen  proportionally  to  their  atmospheric  concentrations.  This  proportionally 
uniform  increase  in  partial  pressures  leads  to  a  proportionally  umform  increase 
in  the  effective  saturation  concentrations,  and  ultimately  a  uniform  increase  in 
the  effective  depth.  Thus,  the  effective  depth  calculated  using  dissolved 
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Table  4 

Results  of  Simultaneous  Methane  and  Oxygen  Transfer  Measurements  at  Various  Hydraulic  Structures  (Uncertainty 
Is  given  to  the  95-percent  confidence  level)  _ 


Table  4  (Continued) 


3i 


(Sheet  2  of  3) 


Table  4  (Concluded) 


WATERrSURFACE  3/2/90 
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Figure  11.  Effective  depth  versus  specific  discharge  at  Rum  River  Dam 
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SPECIFIC  FLOWRATE  (mVs) 


nitrogen  measurements  should  equal  the  effective  depth  calculated  using  dis¬ 
solved  oxygen  measurements. 

Table  5  shows  the  results  of  simultaneous  methane,  oxygen,  and  nitrogen 
transfer  at  the  Rum  River  structure  on  31  August  1993.  At  each  specific  dis¬ 
charge  tested,  the  effective  depth  computed  using  dissolved  nitrogen  concentra¬ 
tions  was  found  to  be  similar  to  the  effective  depth  computed  using  dissolved 
oxygen  data  within  the  confidence  interval  of  the  measurement 
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7  Conclusions 


The  following  conclusions  were  drawn  from  this  study. 

a.  Measurement  of  in  situ  methane  concentration  is  a  viable  technique  for 
determining  the  transfer  efficiency  across  a  hydraulic  structure  or  any 
device  having  a  large  rate  of  gas  transfer  and  a  relatively  short  residence 
time.  It  may  also  have  other  applications,  such  as  measuring  transfer 
efficiency  in  pneumatic  diffusers. 

b.  Methane  is  usually  not  stratified  in  shallow  reservoirs  (less  than  8  m 
deep)  except  under  ice  cover,  where  methane  was  found  to  be  strongly 
stratified. 

c.  Because  methane  stratification  in  deeper  reservoirs  increases  the  uncer¬ 
tainty  in  deteimining  the  mean  methane  concentration  at  the  top  of  the 
hydraulic  structure,  the  methane  tracer  technique  is  most  applicable  to 
low  head  structures. 

d.  Oxygen  and  methane  transfer  efficiencies  were  found  to  be  different 
when  atmospheric  pressure  was  used  to  compute  oxygen  saturation  con¬ 
centration  and  significant  entrainment  of  air  into  the  tail  water  pool 
occurred.  This  difference  is  believed  to  occur  because  entrained  air 
bubbles  pulled  to  a  depth  of  higher  pressure  in  the  taU  water  results  in  a 
larger  saturation  concentration  for  oxygen,  while  methane  saturation  is 
essentially  zero  throughout  This  is  a  source  of  error  in  using  the  trans¬ 
fer  efficiency  of  tracers  to  compute  oxygen  transfer  efficiency  and  in 
using  the  measured  oxygen  transfer  efficiencies  to  compute  the  transfer 
efficiencies  of  other  compounds  of  interest,  such  as  the  many  volatile 
toxic  compounds  in  surface  waters. 

e  'Simultaneous  oxygen-methane  transfer  measurements  and  nitrogen- 
methane  transfer  measurements  have  been  used  to  compute  an  effective 
saturation  concentration  that  is  a  weighted  mean  of  the  samration  con¬ 
centration  that  a  submerged  bubble  experiences  in  its  path  through  a 
hydraulic  structure  and  plunge  pool.  The  measurements  have  also  b^n 
used  to  compute  a  similarly  weighted  mean  depth  of  bubble  penetration, 
called  the  effective  depth,  at  two-gated  control  stractures. 
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/  Generally,  as  discharges  increase,  the  effective  depth  increases,  due 
perhaps  to  the  increase  momentum  pulling  the  bubbles  deeper  in  the 
plunge  pool. 

g.  Measurements  were  used  to  compute  the  effective  depth  as  a  function  of 
specific  discharge  on  three  different  days  at  one  stracture.  At  compara¬ 
ble  discharges,  the  effective  depths  were  found  to  be  similar  within  the 
confidence  interval  of  the  measurements. 

h.  The  effective  depth  computed  from  oxygen-methane  measurements  was 
similar  to  that  computed  from  nitrogen-methane  measurements  within 
the  confidence  interval  of  the  measurements  for  a  range  of  specific  dis¬ 
charges  at  a  stmcture,  thus  verifying  the  concept  of  effective  depth. 

I.  More  of  these  measurements  need  to  be  performed  at  other  stmctures  to 
relate  the  effective  depth  to  parameters  of  the  jet  and  stracture.  Only 
then  can  full  use  of  gas  transfer  measurements  be  made  at  hydraulic 
structures. 
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Photo  2.  Water  sampler 


Photo  3.  Creation  of  headspace  in  sample  bottle 


Appendix  A 

Air  Entrainment  Data 


measured  at  the  Rum  River  Dam  2/17/90 

Qwater  =  108  Cftn/ft 


Distance  from 
plunge  point 
ft 

Run  No.  1 

q.lr 

cfm/ft 

Run  No.  2 

q.lr 

cfm/ft 

0.75 

6.47 

9.53 

2.25 

7.27 

5.73 

3.75 

3.40 

4.20 

5.25 

8.00 

6.47 

6.75 

2.67 

4.20 

8.25 

4.93 

4.93 

9.75 

7.27 

6.47 

11.25 

4.20 

4.93 

12.75 

2.67 

3.40 

14.25 

2.00 

4.20 

15.75 

1.40 

2.00 

17.25 

0.80 

0.80 

18.75 

1 

0 

0 

Total 

76.62 

85.29 

Appendix  A  Air  Entrainment  Data 


Appendix  B 

Description  of  Structures 


Kost  Dam 

The  Kost  Dam  is  located  on  the  Sunrise  River  in  Oiisago  County,  Minne¬ 
sota.  The  dam  consists  of  an  imgated  ogee  spillway.  A  plan  view  and  a  pro¬ 
file  view  of  the  dam  are  shown  in  Figures  B1  and  B2,  respectively. 


Rum  River  Dam 

The  Rum  River  Dam  is  located  on  the  Rum  River,  approximately  500  ft 
upstream  of  the  Main  Street  Bridge  in  Anoka,  MN.  The  structure  consists  of  a 
236-ft-long  fixed  weir,  and  a  20-ft-wide  tainter  gate.  The  tainter  gate  spillway 
is  separated  from  the  main  river  by  a  concrete  pier.  The  plan  view  of  the 
structure  and  the  specific  sampling  locations  are  shown  in  Figure  B3.  The 
profile  of  the  fixed  weir  is  shown  in  Figure  B4.  The  profile  of  the  tainter  gate 
and  spillway  is  shown  in  Figure  B5. 


St.  Cloud  Dam 

The  St.  Qoud  Dam  is  located  on  the  Mississippi  River  approximately 
500-ft  downstream  of  the  10th  Street  Bridge  in  St.  Qoud,  MN.  The  structure 
consists  of  a  500-ft-long  fixed  weir,  and  two  tainter  gates.  The  tainter  gate 
spillway  is  separated  from  the  main  river  by  a  concrete  pier.  The  plan  view  of 
the  structure  and  the  sampling  locations  are  shown  in  Figure  B6.  Only  the 
tainter  gate  portion  of  the  structure  was  investigated  in  this  study.  The  profile 
of  die  tainter  gate  and  spillway  is  shown  in  Figure  B7. 


Smithland  Lock  and  Dam 

The  Smithland  Lock  and  Dam  is  located  on  the  Ohio  River  near  Paducah, 
KY.  A  profile  of  the  dam  is  shown  in  Figure  B8. 
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Figure  B1 .  Plan  view  of  Kost  Dam 


B2 
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Figure  B3.  Plan  view  of  Rum  River  Dam 


Opekiska  Lock  and  Dam 

The  Opekiska  Lock  and  Dam  is  located  on  the  Monogahelia  River  near 
Morgantown,  WV.  A  profile  of  the  dam  is  shown  in  Figure  B9. 
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Figure  B5.  Profile  of  Rum  River  fainter  gate 


FIXED  WEIR 
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Figure  B8.  Profile  view  of  Smithland  Lock  and  Dam 


Figure  B9.  Profile  view  of  Opekiska  Lock  and  Dam 
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Appendix  C 
Notation 


A  Air-water  surface  area 

C  Concentration  in  water 

Ca  Concentration  in  the  air 

Cj,  Time  averaged  methane  concentration  in  a  bubble 

Concentration  downstream  of  the  structure 
Concentration  of  methane  in  a  bubble  as  it  is 
released  to  the  atmosphere 
Equilibrium  concentration  of  methane  in  a  bubble 
downstream  of  a  structure 

Dissolved  nitrogen  concentration 
Saturation  concentration 
C„  Effective  saturation  concentration 

C„  Concentration  upstream  of  the  stmcture 

comp  Molar  fraction  of  gas  in  the  atmosphere 

D  Diffusivity  of  the  measured  compound  in  the  water 

Effective  bubble  depth 
Dq  Diffusivity  of  oxygen  in  water 

E  Gas  transfer  efficiency 

E„  Methane  gas  transfer  efficiency 

Eio  Equivalent  transfer  efficiency  of  oxygen 

F  Water  surface  flux  of  a  volatile  chemical  per  unit  surface 

area 

f-  Indexing  parameter  for  chemical  and  temperature 

fg  Indexing  parameter  for  chemical 

/j  Indexing  parameter  for  temperature 

H  Henry’s  Law  constant 

Liquid  film  coefficient 

Kp  Conversion  constant  from  atmospheres  to  depth  of  water 

M  Molecular  weight  of  methane 

Barometric  pressure 
^H20  Vapor  pressure  of  water 

Partial  pressure  of  nitrogen 
Pq^  Partial  pressure  of  oxygen 

Specific  air  entrainment  rate 
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Specific  flow  rate  of  water 
R  Universal  gas  constant 

t  Residence  time 

T  Water  temperature 

V  Control  volume  over  which  A  and  C  are  measured 

Volume  of  headspace  in  a  sample  bottle 
Volume  of  water  in  a  sample  bottle 
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